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ARBITRARILY CURVED AND TWISTED SPACE BEAMS 

by 

William F. Hunter 
(ABSTRACT) 

A derivation of the equations which govern the deformation of an 
arbitrarily curved and twisted space beam is presented. These equations 
differ from those of the classical theory in that extensional effects 
are included. Other departures from the previous theory are that the 
strain - displacement relations are derived and that the expressions for 
the stress resultants are developed from the strain - displacement 
relations instead of assuming that the resultants are proportional to 
changes in the curvatures. It is shown that the torsional stress resul- 
tant obtained by the classical approach is basically incorrect except 
when the cross-section is circular. 

Using a vector approach the exact expressions for the curvature 
components of a deformed space beam are developed. Because inextension 
of the beam is not assumed an additional term appears in each of the 
linearized curvature expressions. These expressions are utilized in 
the derivation of the strain - displacement relations. The normal and 
shearing physical components of the strain tensor are given. These 
relations are not restricted to beams whose cross-sectional dimensions 
are very small compared to the radius of curvature. Next, a develop- 
ment of the stress resultants is presented. Effects arising from the 


initial twist of the beam are obtained which are not reflected in the 
classical theory. Finally, the six equilibrium equations are derived 
using a vector approach. 

The governing equations are given in the form of twelve first- 
order differential equations. A numerical algorithm is given for 
obtaining the natural vibration characteristics and example problems 
are presented. 
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V. mTRODUCTIQN 

The work presented herein arose from the need to better understand 
the elastic and dynamic behavior of twisted and curved propeller blades* 
It has long been recognized that the stresses, in propeller blades can be 
reduced by designing the blades to be slightly curved so that the centri- 
fugal forces are used to counteract the bending moments arising from 
lift and drag. This technique can be especially beneficial to the design, 
of wind tunnel blades* However, an adequate theory for accomplishing 
this has not yet been developed. Also, there was the need to study the 
effects of curvature on. the vibration characteristics of rotating 
blades. In addition, there is much current interest in analyzing the " 
Oarrieus windmill (sometimes called the vertical -axis windmill ). 

Although these specific problems are not treated in this paper, much of 
the needed analytical development work is presented. . . 

A review of the classical theory for curved and twisted space beams 
as presented by A. E. H<- Love in reference 1 led to an interest to ex- 
tend the theory to include additional effects and to present a fresh 
and more rigorous derivation of the governing differential equations. 

The classical theory was developed by Kirchpff, Cl ebsch, Michel 1 , 
and Basset In the latter part of the last century (in reference 1 Love 
discusses the development and lists references). However, Love is 
responsible for presenting the theory in an organized fashion and his 
work is most often quoted. 


1 
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In the classical theory, the central line. of the beam is assumed 
to be inextensible and the cross-sectional dimensions are assumed to be 
very small compared to the radius of curvature. Love gives a develop- 
ment for the curvature expressions of the central line of a deformed 
beam. Love's work contains an irregularity in that his direction 
cosines do not satisfy the orthogonality relations. The inextensibility 
assumption and the direction cosine anomaly, which are related, have 
been questioned by Waltking (see reference 2) and by Philipson (see 
reference 3). 

In this paper the exact expressions for the curvature components 
of a deformed beam, which is initially curved in space in any arbitrary 
manner, are derived using a vector approach. No assumptions are made in 
the development and, thus, the resulting equations are applicable to 
large, as well as small, deformations. When the three curvature rela- 
tions are linearized an additional term, that Love does not obtain, 
appears in each one. It is shown that these terms arise because exten- 
sional deformation is included. 

A development is presented for the strain-displacement relations of 
space beams. Since the strain distribution is not linear over the 
cross-section, the relations can be applied to beams having much larger 
curvatures than that permitted by the classical theory. The strain- 
displacement relations are used to determine the stress resultants. 

In the classical theory the strain-displacement relations are not 
derived. The stress resultants are assumed to be equal to changes in 
the curvatures times an elastic cross-sectional constant. It is shown 
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that this approach yields an incorrect stress resultant which can lead 
to considerable error. 

A vector derivation is presented for the six beam equilibrium equa- 
tions. These equations are essentially the same as Love's equations 
when linearized for natural vibration solutions. 

The analysis developed for curved beams is applied to natural vibra 
tion examples. The governing equations consist of twelve first-order 
differential equations and the solutions are obtained by using a 
transfer matrix method. Included in the equations are rotary inertia 
and elastic foundation effects. 

Reference 4 is a survey paper by Royster on the subject of curved 
beam vibrations. Practically all studies have been on rings or beams 
whose central lines lie in a plane. In fact, in the sunmary of his 
report Royster states "Only two papers on the vibrations of a curved 
beam of double curvature have been published," Also, regarding solu- 
tions he lists the conclusion "In general, solutions for the coupled 
In plane and out of plane vibration problem, as is the case if the 
central-line of the unstressed beam is not a plane curve* are extremely 
rare." In references 5 and 6 Vol terra applies the Method of Internal 
Constraints to beams of double curvature. He obtains nine coupled 
second-order differential equations having variable coefficients with 
the derivatives being of the nine unknown displacement functions. A 
solution method is not presented. 

In papers pertaining to the vibrations and buckling of circular 
rings the governing equations are arrived at by other methods, as well 
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as from Love’s work, when the deformation is in the plane of the ring. 
However, for coupled torsion and out-of-plane bending deformation the 
equations of Love's classical theory are Invariably used. For such 
deformation the classical theory should be modified because of the pre- 
viously mentioned error associated with one of the stress resultants. 

It is the intent of this dissertation to present a new development 
of the equations which govern the deformation of space beams rather than 
to present a study of the vibration characteristics of curved beams. 
However, a numerical solution method is given for such beams. Compre- 
hensive discussions are presented on the developments of the curvature 
relations, the strain-displacement relations, and the stress resultants 
in view of the classical theory according to Love and other literature. 



VI. LIST OF SYMBOLS 


A 

a, b, c 
[A^] 

[a] 

[B] 
b 

[b] 

S’ 

[C] 
d. d‘ 

CD] 

“^x’ '‘y’ 
dM 

dP 

dS, dS' 


d$, ds* 


Cross-sectional area of beam 
Quantities defined by equations (28) 

Transfer matrix defined by equation (145) 

Matrix appearing in equation (132) 

Cross-sectional constant defined by equations (115), 

(i = 1, 2, ... 9) 

Product of the transfer matrices as given by equation 
(147) 

Bi normal vector of the space curve formed by the axis of 
undeformed beam 

Matrix appearing in equation (132) 

Elastic foundation constants for linear displacements 
Matrix defined by equation (148) 

Modified Darboux vector for undeformed and deformed 
beam, respectively 
Matrix defined by equation (149) 

Elastic foundation constants for rotational displacements 
Mass of differential lengths ds and ds' of beam 
Force acting on differential length ds* 

Lengths of an arbitrary Incremental line element before 
and after deformation (see equation (11)) 

Differential lengths along elastic axis before and after 
deformation 
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i fo » n> 
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d§, d§‘ 
E 




F 





[H] 

^xx* ^yy 


1 

[I] 

J 



Differential lengths defined by equations (80) and (81) 
Elastic modulus 

Normal components of the strain tensor 
Shearing components of the strain tensor 
Orthonormal vectors at point on elastic axis of undeformed 
beam 

Orthonormal vectors at point on elastic axis of 
deformed beam (defined by the rotation of the imbedded . 
!y* vectors) 

Vector of internal forces 
Components of F 
Shear modulus 

Metric coefficients for deformed and undeformed beam, 
respectively 

Matrix defined by equation (140) 

Moments of inertia about the cross-sectional axes 
Polar moment of inertia of cross-section about elastic 
axis 

Denotes station along beam 
Identity matrix 

Saint-Venant torsional stiffness constant defined by one 
of equations (115) 

Components of curvature of the undeformed beam 
Components of curvature of the deformed beam 
Direction cosines defined by reference 1 
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R, R' 


Direction cosines relating e^ to e^, ®y» 

(i = 1,2, 3) 

Direction cosine of the angle between the z coordinate 
line and e^ 

Vector of the internal moments 

Components of M 

Mass per unit length of the undeformed and deformed beam, 
respectively 

Normal vector of the space curve formed by the elastic 
axis of undeformed beam 

Denotes position of point before and after deformation, ; 
respectively . 

Applied distributed force vector per unit length of 
deformed beam 

Components of p 

Applied distributed forces per unit length of undeformed 
beam 

Applied distributed moment per unit length of deformed 
beam 

Components of q 

Applied distributed moments per unit length of undeformed 
beam 

Radius of curvature 

Position vectors to point on elastic axis before and after 
deformation, respectively 
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r, r 


s, s' 


t 


[U] 


U, V, w 


iJ 


X. y 

{Y} 

{Y'} 


z 



Y 

AR 

As' 

e 


Position vectors to point in cross-section before and 
after deformation, respectively 
Coordinate along elastic axis of undeformed and deformed 
beam, respectively 

Unit tangent vector of the space curve formed by the 
elastic axis of the undeformed beam 
Matrix defined by equation (151) 

Displacements of the elastic axis in the e , e , and e 

-•X 'rj ^ 

directions 

Second partial derivative of u with respect to time 
Saint-Venant warping function 
Cross-sectional coordinates. 

Column vector of the variables of the vibration problem, 
defined by equation (138) 

Column vectors of the variable derivatives, defined by 
equation (139) 

Coordinate line defined by the intersection of the 
coordinate surfaces x = constant and y = constant 
Quantities defined by equations (20) 

Direction cosine of reference 1 

Orientation of e^ relative to n (see figure 2) 

Strain tensor 

Displacement vector of a point on the elastic axis 

Incremental length of deformed beam 

Tensorial extensional strain of the elastic axis 



Eulerian rotation 
Curvilinear coordinate 

Angle between the and coordinate lines before 
and after deformation, respectively 
Curvature of the elastic axis of the undeformed beam 
Mass density of beam 

Normal stress in direction of the 2 coordinate line 

Torsion of the undeformed beam 

Rotational displacement of the beam about e^ 

Eulerian rotation 
Natural vibration frequency 



VII. ANALYSIS 


In this chapter developments will be presented for curvature 
relations, the strain-displacement relations, the stress resultant 
expressions, and the equilibrium equations for arbitrarily curved and 
twisted space beams. Each of the first three developments is followed 
by a discussion of the development.. These discussions are included in 
this chapter since they are an in-depth analysis of the developments 
in relation to previous work instead of being of a general discussion 
nature. Besides being convenient to the reader, it is thought that 
the placement of the discussions in this chapter will give the reader^ | 
a better understanding of the developments which follow each. 

Also, a summary of the governing differential equations for natural 
vibration is given. This is followed by a description of a numerical 
method for solving them, , 

* » ' f - 

A. Development of Curvature Relations. 

Consider a curved beam whose elastic axis forms any general space 
curve. The elastic axis of the^ beam before and after deformation is 
shown in figure 1 where primes denote the deformed state: In this 

figure, R is the position vector to point ' P on the elastic axis, 
s is the coordinate along the elastic axis.i . AR is the displacement 
vector of point P, and e^^, e^, and e^ are orthonormal vectors with 
e^ tangent to the curve s. 

The unit vectors e^ and are taken to be aligned with the 
principal axes of the beam's cross-section. The orientation of e 
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and 6y Is specified relative to the space curve's normal, n, as shown 

in figure 2. The vectors t, n, and b are the tangential , normal, 

and bi normal unit vectors, respectively, of the space curve. They 

are related to e , e , and e by the transformation 
-X •^y ~z 

n = e cos Y - e sin y 
^ ~x ~,y 

‘ b = e sin y + e cos y 


The shape of a space curve can be completely described by two 
parameters: the curvature k(s) and the torsion t(s). The deriva- 

tives, of t, n, and b with respect to s are given by the well 
known Frenet-Serrett formulas as 



Similar expressions can be obtained for the e , e , e_ trihedron. 

~x ~y -Z 

They are most simply obtained by modifying the Darboux vector {see 
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reference 7 or any other book on differential geometry) and operating 
with it. The Darboux vector is the rotation rate of the trihedron at 
point P as P moves along the curve at unit velocity. Or, it can be 
viewed as rotation per unit length of the curve. Since y is not taken 
to be constant, the modified Darboux vector is ' < .. 

d = (t + ^)t + Kb > (3) 

Substituting equations (1) into the above gives 

^ ~ + K sin Y e^ + K cos y e (4) 

~ ds ~z ~x 

The quantities (t + ^), (k sin y), and (k cos y) are known as the 
components of curvature. Letting 

k = K sin Y 

ky = K cos y 

k “ T + ^ 
z ds 

the modified Darboux vector may be rewritten as 
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In the manner used in dynamics to obtain the derivatives of unit 
vectors with respect to time, it can be shown that the derivatives of 
e^> e^, and e^ with respect to s are 




Hence , 




It is noted that in the preceding equations, k , k , and . k are 

X y i 

the curvature components before deformation. The goal is to arrive at 
expressions for the curvature components of the deformed beam in terms 
of the elastic displacements. 

As shown in figure 1, position vector to the displaced point P' 
is ' 
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■ R' = R + AR ■ ■ (9) 

where 

AR = u e + V e + w e ’ ; (10) 

~ -A -J ~Z ’ 

and u, V, and w are the elastic displacements of point P in the 
directions of the reference unit vectors of the undeformed beam. 

If dS is the length of an arbitrary incremental line element in 
an elastic solid and dS' is its length after deformation, then (see 
reference 8 or other book on elasticity or continuum mechanics) 

(dS')^ - (dS)^ = 2 Y^jde’de^ (i, 0=1, 2. 3) (11) 

where the y.. are the components of the strain tensor and the o’ are 

IJ 

the curvilinear coordinates. 

Applying this equation to a curved beam with the . incremental line 

element being chosen to lie along the elastic axis such that 
1 ? 3 

de = de = 0, de = ds, and Y33 = ^ gives 

(ds' - (ds)^ = 2 e ds^ (12) 


Or, 
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„ 1 r ' \ T -| 


(13) 


where e is tensor ial extensional strain of the centroidal axis. Since 


(ds‘ = dR' ■ dR' 


(14) 


the above may be rewritten as 


dR‘ dR‘ 


c = T [ 


2 '-ds ds 


- 1 ] 


(15) 


Differentiating equation (9) with respect to s. 


dR' dR , d(AR) 

ds ds ds~ 


06 ) 


Also, it is noted that 

dR 

Substituting the above into equation (16) and the result into equation 
(15) leads to 


d(AR) , d(AR) d(Aj^) 

® ds !z 2 ds ds 


(18) 
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Differentiating equation (10) with respect to s and applying 
equations (8) yields 

d(Aj?) 

~di 





which is the nonlinear expression for the extensional strain of the 
elastic axis. The significance of the“ above result in view of 
previous work will be discussed later. 

Equation (16) may be rewritten by substituting equations (17) 
and (19), 


dR’ 

dT" "x!x ^ “y!y ^ ^ (22) 
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From equation (12) it is noted that 


ds‘ = (1 + 2e)^^^ds 


(23) 


Hence, 


dR' . dR' 

g^= (1 - 


(24) 


Also, 


dR' 

^ " 5z 


(25) 


Combining equations (22), (24), and (25) results in 


e^ = (1 + 2E)-l/2|-^^g^ ^ ^ ^ 


x^x y^y 


z''2- 


(26) 


The above equation gives the direction of the unit vector tangent 
to the deformed elastic axis in terms of the displacements, initial 
curvature, and unit vectors of the undeformed rod. 

Equation (26) may be differentiated with respect to s' by 
applying equation (23) and substituting equations (8). Doing this 
yields 


de ' 

"r = (1 + 2e)"^ (a e^ + b e^ +.c e^) 


ds 


(27) 



where 



a = -ayk^ + (1 + ^ ^ . 0 + 2,)-1 a. ^ 

da -1 j 

b = - (1 + a^)k^ + - (1 + 2e)‘‘ ay g| 

c = -a^ky + «yk^ + dT - (1 + 2 e)'^ 0 + a^) ^ 

de 

and from equation (21), ^ is given by 


dg _ 
ds 



da 

“y V“z> 




The direction cosines which define the orientation of the trihedrons 
of the undeformed and deformed beam relative to each other are identi- 
fied in the table below. 



the displacements by equation (26) and are 
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1, = (1 + 2e)-’''2 a. 


nij = (1 + 2s) a. 


n, = (1 + 2s)'^''^(V + a ) . 




)( 30 ) 


In the same manner as before, a modified Darboux vector may also 
be written for the deformed beam as 


d' = k'e' + k'e' + k‘e‘ 
^ xSx y^y z~z 


(31) 


Using this vector, the derivatives of the unit vectors of the deformed 
beam are found to be 


de' 
de! 


A 


J. = -k'e* + k*e* 
ds' Slyi \?2 


de* 

= k'e' - k'e' 
ds’ ’^y!x *^x!y 


\ (32) 


J 


From the third of the above equations and equation (27), 


kyO^ “ Kty " ^ 5x ^ ?y ^ 5z^ 
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The transformation between the unit vectors in terms of the direction 
cosines is 


?x = ^ ^2!y " ^3 !z 


■!y = "’r!x ^ "’zS; ^ "'sS; 


Sz = "l!;; ^>'2!; ^ "3 !z 


) (34) 


Substituting the above into equation (33) and rearranging gives 


-1 




+ (al2 f cn2)ey 


+ (al 2 + bm^ + cn 2 )e^] (35) 


Thus, 


~ ^ ^(^12 ’’’ ^'^2 ^*^2^ 

ky = (1 + 2e)”^{al^ + bm-j +.cn^) 

0 = (1 + 2 e )~^( al 2 ^ ^*^3 ^^ 3 ^ 



•) ( 36 ) 

j 
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Two of the curvature components have now been expressed in terms of ,, 
a, b, and c (defined by equations (28)) and the direction cosines,. 

It can be shown that the third of equations (36.) is exactly satisfied 
by the previously defined a, b, c, and- n^. Also, by squaring 

equations (36), adding, and applying the transformation orthogonality, 
relations it is seen that 

+ k'^ = (1 + 2E)'^(a^ + + ;c^) , (37) 

X . y 


From the first of eqs. (32), 




1 

(38) 


Using equation (23) and substituting the transformations 




= ^ "’l!y ^ "l!z 


^ ^ ^ "’25y ^ "2?z 


> (39) 


into equation (38) gives 


k' = (1 + 2e)-^''^[^ (l^e^ + m^e + n,e^)] * [1 + m2!. +' >26^] 


Differentiating the above and substituting equations (8), for the 
derivatives of the unit vectors results in 
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dl , dm, dn, 

+ (nijl^ - l2'"l)'^z ^2 3s“ '"2 3i~ * "2 dF"^ 


(40) 


The table of direction cosines given on page 20 is an orthogonal 
matrix. The inverse of an orthogonal matrix is equal to its transpose 
and its determinant is equal to unity. Thus, each element of an 
orthogonal matrix is equal to its cofactor. Applying’ this theorem to 
the last row of the direction cosine matrix yields the following 
orthogonality relations 




1 2 m^ ^2 ^ ^2 ^1 


*^1 *^2 ^2 '"l 


7 


03 1^ m^ - I2 ni^ 


The above relations may have been obtained by substituting the trans- 
formation expressions for e' , e', and e' into e‘ x e' = e' and 

~x •'.y ~z. --X' ~y ~z 

carrying out the operation as Novozhilov did in reference 9. Equation 
(40) can now be rewritten as 


” (1 + 2e) + m3k^ + n3k^ + I2 


+ m 


dm: 
2 dT 


dn 


+ r>2 5^) (42) 



25 


Equations (36) and (42) give the curvature components in terms of 

the direction cosines, the initial curvatures, and the quantities a, 

b, and c. The direction cosines m^, and n^ have been expressed 

in terms of the displacements u, v, and w. These cosines (1^, 

and n^) define the orientation of e^ relative to the trihedron of the 

undeformed beam (i.e., e , e , and e ). Besides knowing the orientation 

~x ~y ~z 

of e^, one additional angle (or direction cosine) is needed to fully 
specify the orientation of the trihedron of the deformed beam relative 
to the trihedron of the undeformed beam (in other words, the orientation 
of e^ and needs to be specified). This additional rotation 
joins u, v, and w as one of the elastic displacements (or variables) 
of the problem. The direction cosines ni^, n^, 12* m^, and n2 
can be expressed in terms of the known cosines (l^, and n^) and 
the new elastic displacement which has not yet been -defined. 

When the beam is deformed, the principal flexure-torsion trihedron 
at any general point on the elastic axis undergoes translation and 
rotation. The rotation may be expressed in terms of Eulerian angles. 

It has been found that the most suitable Eulerian angle system for this 
application is that used in aeronautical and aerospace engineering 
(see reference 10). For this system all three Euler angles are small 
if small deformations are assumed. This is not the case for the Euler 
angle system of reference 11. 

Let xyz be the initial reference frame aligned with the e^, e^, 

e^ trihedron and let x’y'z' be the final reference frame which is 

aligned with the e‘, e', e' trihedron. The three rotations will be 

~x ~y ..z 
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taken in the following order:. 

1. A positive rotation ij; about y axis resulting in XYZ reference 

2. A positive rotation 0 about X axis resulting in X'Y'ZV 

reference 

3. A positive rotation (j) about V axis resulting in x‘y'z' 

reference 

It is noted that the above is different from reference 10 in that the 
axes are not named the same. The three transformations are 



r'' 







f ^ 


X 


COS 


0 

-sin 


X 


Y 

> = 

0 


1 

0 

< 

\ 

y 


z 


sin 


.0 

cos Ip 


2 



/ 






^ J 

c 

'N 





— 


r 'i 

1 


> 

1 


0 ■’ 



X 

' 


0 


cos e 

sin 8 






0 


-sin 8, 

cos 6 

1 

z 

V. 

J 







^ J 

r 



- 





r ^ 




COS 

4) 

sin $ 

0 


X' 




-sin 


cos‘ 4 

P 

< 

Y' 


z' 


0 


0 

1 


V 








_ 

< J 


Substituting and performing the matrix operations yields 
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y 

"iT 

z' 


V J 



cos (p cos Ip sin <p cos 0 -cos 4> sin \p 

+s1n (J) sin B sin i/;* +sin sin 9 cos ip 

-sin (p cos ^ cos (p cos 0 sin (p sin 

+COS (p sin 0 sin ip +cos (p sin 0 cos 


cos 0 sin 


-sin 6 


cos 6 cos 




I y \ (43) 


z 

V. y 


The above gives the previously discussed direction cosine matrix 
in terms of Q, and <|). The angle <p will be taken as the needed ad- 
ditional displacement variable. The angle (p is the rotation about the 
z' axis (or elastic axis) that brings the X' and Y' axes to the 
final x' and y‘ position. 

From the direction cosine matrix of equation (43), 

1^ = cos 0 sin yp 


= -sin 0' 


Or. 


n^ = cos 0 cos i|) 


0 = -sin~^ m. 


-1 ^3 
= tan ‘ (^) 

"3 




> ( 44 ) 



.1 
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By substituting equations (30) the above expressions for 0 and ^ may 
be rewritten in terms of the displacements as ' 



Using the direction cosines defined by equation (43), the expressions 
for and given by equation (36) become 


1 2 ^ [(-sin (p cos + cos i|) sin 6 sin 4>)a + (cos cos 6)b 

+ (sin (p sin + cos 4) sin 6 cos i{i)c] 

1 

' 2 ~ £ [(cos 0 cos ip + sin (p sin 6 sin 4;)a + (sin > cos 6)b 

+ (-cos <p sin + sin 4> sin 8 cos i|;)c] 
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Substituting the above into equations (46) gives 

= (1 + + 2e - + 2e)^^^(l + a^)sin (p + cos (j>]a 

- [(1 + 2e - ay)cos <p]b + [-(1 + 2e)^^^a^ sin <j) + ay(l + a^)cos c))]cj- 

(48) 

= (1 + 2e)'^/^(l + 2e - + 2e)’^^(l + a^)cos # - sin >]a 

+ [(1 + 2c - a^)sin <l)]b + [-(1 + 2e)^^^a cos iji - a (1 + oi )s1n <(i]cf 

(49) 

The expression for in terms of the direction cosines is given 
by equation (42). The last three terms of this equation may be expressed 
in terms of the Euler angles by substituting the direction cosines as' 
defined by the matrix of equation (43). Performing the operations and 
combining terms yields the surprisingly simple result 


1 ^ m ^ + n ^ ^ - sin e ^ 
h ds * "’2 ds ^ "2 ds ds ® ds 


(50) 


Substituting equations (30) and (50) into equation (42) gives 


= (1 + 2e) + k^a^ + k^(l + a^)] + (1 + 2e)“^^^.(^ - sin 6 ^) 

. (51 ) 

The last term of equation (51 ) may be rewritten in terms of the dis- 
placements as 
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m. 


dl 


dn. 




2 . ,2 ^"3 ds 3 ds 
"3 ^ ^3 


Using equations (30) it can be shown that 


(52) 


dl~ dn^ , da da 

"3 ^ ■ ^3 dT ° * 2e)-'[(l + aj - a, (53) 


z’ ds X ds 


Substituting the above and equations (30) into equation (52) gives 


•sin 8 ^ = 


(1 + 2e)'^''^a 


(1 + 


4[(1 + «J 


da 


da. 


- a. 


1 * ds X ds 


(54) 


Noting from equation (21) that 


(1 + a^)^ + =1 + 2e - a^ 


and substituting equation (54) into equation (51) gives the third 
curvature component as 


a da > 


1 + 2e - a 


(55) 


Equations (48), (49), and (55) are the equations for the curvature 
components of the elastic axis of the deformed beam in terms of the 


31 


elastic displacements and the initial curvature. The quantities 

“x’ “y’ a. b, and c appearing in the equations are defined as 

functions of the displacements u. v, and w by equations (20), (21), 
and (28). 

It is noted that the expressions derived for k', k', and k' are 

X y z 

exact. No assumptions whatever have been made. 

Also, it is observed that different forms of the curvature compo- 
nents can be obtained by altering the order of the Eulerian rotations. 
For example, if the rotations are taken in the order 

1. A positive rotation 0 about x axis resulting in XYZ reference 

2. A positive rotation ip about Y axis resulting in X'Y'Z' > 

reference , , 

3., A position rotation about Z‘ axis resulting in x'y’z’ 
reference 

then the curvature components are given by 
K; = (1 + 2 c)- 3/2 (i + 2 , - a2 )-l/2 j^^, ^ 

+ [-(1 + 2e)’^^(1 + a^ycos 4, - sin <{,]b + [{1 + 26)^''^^^ cos 4 

- 0^(1 + a^)sin ^]c} (5gj 

k; = (1 + 2 c)- 3 / 2 (i + 2, . a2)-V2{[(i + 3^ . 

+ [(1 + 2e)’''^(1 + a^)sin cos (f.]b + [-(1 + 2e)’^2 ^ 

- “jjO + a^) cos (f]cl 


( 57 ) 
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k' = (1 + 2e)'''{k^a^ + kyOy + k^O + c.^) + 0 + 2e) 


1/2 ^ 
ds 


a. 


da aa X 

^ “z^ 55^ ■ “y rf 


da. 


(1 + + Qy 


( 58 ) 


As a check on the expressions given for k' and k‘ it can be 

X y . 

shown that they do satisfy equation (37). 

The exact curvature expressions, which are highly nonlinear, may 
be approximated to any desired degree by applying the binomial theorem 
and substituting the trigonometric expansions for sin <|> and cos 
Using equations (48), (49), and (55), the nonlinear approximations which 
contain terms of order no higher than products of two displacements (or 
squares of displacements) are 


da. 


‘‘x ° '‘x ■ ‘‘z“x ' '^x“z * v ' di^ ■ ''x“x ■ 2 ‘'x"^ * '^x“z ' 2 


da. 


da 


+ 2k,a„a, - k.cx,,* - k„a,iti + a„ 7;;^ + 2a^ + if. 


z x z z y z 


y ds 


da^ 

HF 


(59) 


k'=k “ka “ka - k4) + - .1 k a^ - k a^ + k di^ - i k 6^ 

y y 2 y y z ds 2 y x y y y z 2 y^ 


da da da 

-kaa +2kaa + ka6 + kad) - a ~ 2a -r^ + (t (60) 
x x y z y z z x^ x z^ x ds z ds ^ ds ■ ^ 


k‘®k + ka + ka - ka +^-ka^-ka^ + ka^ 
z 2 X X y y z z ds z x . z y z z 


da j . 

-2kaa -2kaa + a j, ^ - a 
x X z y y z y ds z ds 


( 61 ) 
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It is seen from the above that the linear approximations for the curva- 
ture components are 


da 


V ■ di^ ■ 

(62) 

da 


- "x'*’ ^ ds 

(63) 

ds 

(64) 

the Eulerian rotations 

has no 


effect on the linearized equations so long as there is a rotation about 
each of the body axes. Also, it is noted that the linear relations can 
be obtained most easily from equations (46) and (51) by making small 
angle assumptions for ij;, e, and 4>, 

Using equations (20) the linearized curvature components may be 
expressed explicitly in terms of u, v, w, and 4» as 




p p dk 


k- =k^.^.2k^^.(k2-k2)u 


dk dk^ 

- * di-)'' M 


(65) 


( 66 ) 
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k; = k, + k„ # + k. . ^ - k 


dv 


^ dw 

X ds 'V ds ds 


+ ^ + 2kyk^u - 2k..k_v 


X z 


(67) 


Discuss ton of Curvature Relations 

A development is presented for the exact expressions of the curva- 
ture components of a deformed beam which is Initially curved In any 
arbitrary manner. The exact equations have not been derived heretofore. 
The curvature relations are developed using a vector approach. In con- 
trast to nearly all previous work the beam is not assumed to be inex- 
tensible. The principal cross-sectional axes are assumed to have a 
variable orientation relative to the. curve’s normal and blnormal. Also, 
nonlinear and linear approximations of the curvature equations are given 
In reference 1, Love gives expressions simlTar to the linearized 
equations presented herein. However, there is one term in each of 
equations (62), (63), and (64) that Love does not obtain. The presented 
vector derivation yields these additional terms which are apparently of 
the same order of magnitude as the other terms. The pew terms In the 
k;. k^, and k; expressions are (-k^ a^), .(-k^ a^), and (-k^ a^), 
respectively. 

There are two aspects of Love’s work on curved beams which have 
been questioned previously. First, there is the anomaly in Love's 
work that his direction cosines do not satisfy the orthogonality 
relationships. In particular, the sum of the squares of the direction 
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cosines relating e' to e , e^, and e is not equal to unity. 
Secondly, the validity and applicability of Lovers inextensibility 
relation for curved beams is questioned. 

In the presented development all of the direction cosine ortho- 
gonality conditions are, of course, satisfied exactly. For example, it 
is seen that 1| + = 1 is satisfied exactly by equations (30) 

after equation (21) is substituted. 

The vector approach can be used to explain Love's direction cosine 
anomaly which has not been fully understood before. It shows that 
products of deformations must be included in the development if the 
direction cosines are to be consistent. Also, the reason for the exact! 
formulation yielding additional terms in the linearized equations can be 
shown. 

Love denotes the direction cosines of the angles between e^ and 
each of e^, e^, and e^ by and N^, respectively. Expressed 

in the notation used in this paper. Love finds the direction cosines 
to be 

4 ' “x> 4 = “y* 4 = ’ * “z • 

2 2 

It is obvious that these direction cosines do not satisfy L^ + + 

2 

= 1 , but give 

= 1 + 2a + 

333 zxyz 


= 1 + 2e 
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In arriving at his curvature expressions. Love applies an Inextehsi- 
billty condition which he takes to mean that the elastic axis is un- 
strained. It 1s not thought that such a condition is applicable to 
curved beams. In general, the elastic axis will extend since loading at 
one point along the beam will create at some other point an internal 
force component which is parallel to the elastic axis at that point if 
the beam Is curved. In the presented formulation there was no need for 
imposing an inextensibility relation. 

If a curved beam were assumed to be inextensible, then the 
inextensibility condition is given by setting e of equation (21) equal 
to zero. That is , 

== «z ^ "y ^ ^ (68) 

Love deduced the inextensibility condition tb be 



From equation (68) it is seen that the exact Inextbnsibility assumption - 
actually requires that 

< 0 

If the inextensibility condition, c - 0, were applied to equations 
(48), (49), and (55), then the linearized equations become 
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da 

' “x - ^“z ^ "'y^ - di^ 


da. 




y "y '"z^y * Vz ' ds 


k* = k +ka .+ ka + ka + 
z z X X y y z z ds 


By comparing the above with equations (62) through (64) it is seen that 
imposing the inextensibility condition c = 0 changes the sign of the 
terms containing a^. 

From the definitions of a , a , and a (equations (20)) it appears 

X y z 

that the terms k^ a^, k^ a^, and k^ a^ are of the same order of ^ ' 
magnitude as the other deformation terms in the curvature expressions. 
However, if the beam were assumed to be inextensible, then it is seen 
from equation (68). that would be a higher order term since solving 
equation (68) for gives 


az = -^(a^ + aj) (69) 


Hence, if the beam were assumed inextensible then the terms could 

be neglected and the linearized curvature equations reduce to Love's. - 

But there is no reason to impose the inextensibility condition; the 

curvature equations can be derived without using it. Also, the elastic 

axes of curved beams do, in fact, extend. In addition, a^ is 

apparently of the same order as a and a (see equations (20)). 

X y 

Moreover, the inextensibility assumption imposes a constraint relation 
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between u, v, and w as is given by equation (69). 

Using the notation of this paper. Love gives the direction cosines 


as 


■t- 

V 


■S*v- 

k^w 

It 

•5T-V 

'<x'' 

(= 1 + a^) 


and states 

2 2 2 

"The equation = 1 leads, when we neglect squares and 

products of u, V, w, to the equation 


^ - k u 
ds V 




= 0 


which expresses the condition that the central-line is unextended. In 
consequence of this equation we have = 1“. 

It appears that Love creates his inextensibility condition to 
explain away the anomaly that his sum of the cosine squares is not 
equal to unity. Regardless of the correctness or incorrectness of 
Love's inextensibility relation, the manner in which it is deduced 
does not seem to be a logical method for arriving at it. It should come 
from a strain expression. Love should not have concluded from 
L^ + = 1 that the elastic axis of a curve beam does not 
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extend. The problem was that his direction cosines were in slight 
error. 

It is interesting to note that the direction cosines developed 
herein satisfy ^3 ^ the inextensible case as well as 

for the extensional case. 

It was mentioned earlier that Love's direction cosine problem 
arises because products of deformation are not considered. This can 
be illustrated by using the presented vector approach with Love's 
assumptions applied. Neglecting products of deformation in the 
extensional strain expression of equation (18) gives 

d(AR) 

e = — • e 

ds -2 

d(AR) 

where is given by equation (19). Substituting yields 


The above result is, of course. Love's inextensibility condition when 
set equal ,to zero. Equation (22) is repeated below for convenience. 

dR' 

For inextension of the elastic axis, ds' = ds. Thus, 

dR' dR' 
ds " ds' ?2 



and the above becomes 


e'=ae + ae +(1+a)e 
-z X-S.X y~y ~2 


The above is Love's result In vector notation since the direction 
cosines defined by the above are the same as Love's. Next, Love reasons 
that ~ 0 (inextensibility condition) and the above becomes 


e'=a,e + ae + e 
^2 x^x yZy ^2 


Although the orthogonality condition on the direction cosines is still 
not satisfied exactly, it is now satisfied to a higher order when 
=0 is assumed. It is seen that Love's inextensibility condition 
(a^ = 0) and his direction cosines are a consequence of neglecting 
products of deformations and assuming ds' = ds. It is noted that the 
latter expression above for e^ is the same that Ojalvo and Newman 
obtained in reference 12. Ojalvo and Newman derived Love's expressions 
for the curvature components in the appendix of their paper. 

A close examination of the presented derivation shows that the 
additional linear terms are due to the combination of not setting 
ds' = ds and not setting = 0. Consider, for example, how the 
term {- k a^) arises in the k* relation. In the second of equations 
(46) there appears the product 



41 


where the factor inside the brackets comes from the definition of a. 

The term (1 + 2e)’^ arises from setting ds* = (1 + ds when 

derivatives were taken. The fourth of equations (47) gives the 
expression for cos in terms of the displacements. If the binomial 
theorem is applied and equation (21) is substituted > it is seen that 
cos does not contain a linear term. Assuming <p to be small the 
linear approximation of the above product is 

0 + 2e)'^ cos (f cos i)< [(1 + o^)k ] » (1 - 2e)(l + ct^)ky 

= (1 - 2a^)(l + a^)ky 

= ^ ■ Vz 

Thus, it is seen that this additional linear term, - k^ a^, comes from 
setting ds' = (1 + 2c) ' ds and not letting = 0. If equation 

(63) is considered to be the correct linearized equation for k^, then 
the assumptions of ds' = ds and a = 0 cause an error in k' which 
is equal to (+ k^ a^) since the term does not appear when such assump- 
tions are made. The additional linear terms do not arise from the 
Euler ian rotations. 

Another aspect of the development which should be mentioned con- 
cerns the rotation angle <p. In the present derivation the angle cj) 
is precisely defined as the Eulerian rotation about the body axis 
that is tangent to the elastic axis. Love's equations contain the 
rotational displacement 3 which, evidently, is equivalent to the 
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angle when the exact equations are linearized. Love defines 3 to 
be the direction cosine of the angle between e' and e , From Love’s 

—A 

direction cosines (L-| = 1, M-j = 3» L 2 = -3* M 2 = 1) it is seen that 3 
is essentially a small rotation about the elastic axis with the 
direction of the elastic axis unchanged. Other authors (see, for 
example, references 3 and 12) have followed Love's work in defining 
the rotation angle 3. However, this definition is not satisfactory 
when considering the exact expressions or the nonlinear approximations 
for the curvatures. 

Now consider an untwisted circular beam of radius R whose normal 
is in the e direction such that k - 1/R and k = k = 0. The 
linearized curvature of the deformed beam as given by equation (66) 
reduces to 

k: • i + ^ + V (70) 

The above is the well known formula for the curvature of a circular 
beam deformed in plane (see, for example, reference 13), This result 
is often associated with (or identified as) inextensible deformation. 
However, since the above is not based on inextension the above formula 
should be properly regarded as the linearized extensional curvature 
relation for a deformed circular beam. 

Aside, it is noted that if the beam formed a planar curve instead 
of a circle then the linearized curvature would be 



43 



+ ^ + 
ds2 




”dT 


Love arrives at the same expression (equation (70)) for the 
curvature of a circular beam and this seems surprising at first in 
view of the fact that equation (70) was obtained utilizing the addi- 
tional term which Love's equation for does not contain. The 

reason why Love obtained the same result can be explained. Using the 
notation of this paper. Love's curvature expression for the circular 
beam reduces to 


k' 


1 doc 
1 + 

R ds 


, 1 ^ d^u ^ 1 dw 


ds' 


R ds 


(71) 


Love's inextensibility condition reduces to 


ds R 


or 


dw _ u 
ds “ R 


which he substitutes into k^ when considering applications. This is 
the same as subtracting ^ (^ “ ^y right hand 

side of equation (71). This, of course. Love can do since ® 

from his point of view (but from the extensional point of view this 
could not be done since is not equal to zero and is not of higher 
order). Subtracting (^ - ^), or k^ from equation (71) gives 
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equation (70). It was pointed out earlier that if the extensional 
linearized expression for is considered to be correct then Love's 
expression for would be in error by the amount (+ k^ a^) (this is 
not saying that Love's work is in error - he' assumed inextension). 
Continuing to examine Love's work from an extensional point of view, 
the subtraction of ^ - ^) causes an error of (- k^ a^). Thus, the 

two errors cancel each other and this is the reason the curvature results 
are the same. 

In equations (65) through (67), the linearized curvature components 
are given in terms of u, v, w and 4* Love's inextensional results 
are presented below in the same notation for comparison. 


k'=k - ^-2k ^+k 

•^x "^x ds ’^x ds ds “ 


dk 


9 y 

+ k^v + {-kyk^ + 3^)w + ky4> 


2 

If ' = k + ^ ■ - 7k + |f 

S > ^ ^'^2 ds "'y ds '^z*' 

dk. dk 

- dT ^ di^>« - '‘x* 


= "^z '‘x f * '‘y W ■ 


(72) 


(73) 

(74) 


Depending upon the problem, one applying the first two of Love's 
curvature equations may likely apply Love's inextensional relation 
which is 
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_ dw f. ,, 
“z - d? • 




= 0 


But it is seen that subtracting from equation (72) and 

subtracting k a from equation (73) yields equations (65) and (66), 

respectively, which are the linearized extensional results. Also, 

subtracting k^ from equation (74) gives equation (67)i however, 

one using equation (74) would not be expected to do this since the 

dw 

equation does not contain a ^ term. 

The reason for Love’s expressions for k^ and k^ (after Love’s 
inextensibility condition is applied) being identical to the linear 
extensional expressions is the same as that given earlier for the 
circular beams. For extensional behavior. Love's equation for k^, . 
for example, is in error by the amount k^ a^. Then, when Love's 
inextensibility condition is applied to his equation, k^ is reduced 
by the amount k a , which is not zero from the extensional point of 

y 2 

view. 

t 

It has been shown that Love's expressions for and k^ with 

his inextensibility relation applied are the same as the linearized 

extensional relations for k’ and k* Thus, authors who have used 

X y 

Love's equations (with his inextensibility relation applied) in the 
past have mistakenly assumed that an inextensibility restriction was 
imposed and that their results described inextensional behavior. 

It is noted that although the longitudinal strain e of the 
elastic axis is defined using the definition of the nonlinear strain 
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tensor for curvilinear coordinates (see equation (11)), the knowledge 
that e is a strain is not needed in the development since c is 
eventually expressed in terms of the displacements (see equations (21) 
and (20)). The derivation of the curvature relations is purely a 
geometry problem. The development could have begun with the quantity 
given by equation (13) being denoted as e. Although e does appear 
in the final expressions for the exact curvature components, it may be 
eliminated by substituting equation (21). For the curvature expressions 
to be mathematically correct the extensional strain c must not be set 
equal to zero. 

C. Development of Strain^Displacement Relations 

The strain-displacement relations will be developed with the 
following assumptions being made regarding the displacement of points 
within the cross-section. 

1. Torsional deformation causes warping of the cross-section and 
the out-of-plane displacements are proportional to the Saint-Venant 
warping function. 

2. The cross-sectional body-axis coordinates of points off the 
elastic axis do not change during deformation. This assumption is in 
agreement with Saint-Venant torsion theory that cross-sections do not 
change in size or shape during deformation. 

3. Shear deformation is a negligible effect. Thus, except for 
warping due to torsion, cross-sections remain normal to the elastic 
axis during deformation. 
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Letting x and y be the cross-sectional body-axis coordinates 
of a point in the beam and letting r and r* be position vectors to 
the point before and after deformation, respectively, then 

r(x,y, s) = R(s) + xe^(s) + yey(s) (75) 

r'(x,y.s*) = R'(s') + xe;(s') + ye;{s‘) + W(x,y) ^|^e;(s') 

(76) 

where W(x,y) is the Saint-Venant warping displacement function. Taking 
the differentials of the above relations and substituting equations (8), 
(17), (25), and (32) for the derivatives of the vectors gives. 

dr = (dx - yk^ds)e^ + (dy + xk^ds)e^ 

+ (1 - xky + ykj^)ds e^ (77) 

dr' = (dx - yk^ds' + Wk^ ||t- ds.‘ )e^ 

+ (dy + xk^ds' - Wk^ ds' )Cy . 

+ (ds'-xk;ds'+yk;;ds' 

(78) 

The magnitudes of dr and dr' are the lengths of an arbitrary 
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differential line element in the beam before and after deformation* 
respectively. If the differential length is denoted by d§, then 


(d§)^ = dr • dr 
(d§')^ = dr’ • dr’ 



Substituting equations (77) and (78) into the above and replacing ds'^ 
by 

I 

ds' = (1 + ds , 

I 

in the second of the resulting equations and then collecting terms 
yields 

(d§)^ = dx^ + dy^ + (1 - 2xky + 2yk^ + x^k^ - 2xyk^k^ 

+ + x^k?)ds^ - 2yk dxds + 2xk dyds (SO) 
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(dS')2 = [1 + (|^)^0^{l + 2e)-l]dx^ 

+ [1 + + Ze)’^ V 

+ [1 + 2e + y^k'^n + 2e) - 2yWk^k’ ^ 0 + Ze)^''^ 

+ + x^k'^0 + 2e) - 2xWk'k^ ^ 2e)’''2 

2 2 
+ W^k;^(^) - 2xyk‘(l + 2e) + 2ykj;(l + 2e) + 2W ^ 

^ « ds 

+ xV^(l + 2e) - 2xykik‘0 ,+ 2e) - 2xWk‘ — y 
V ^ ds 

2 2 ^ 

+ y^ky (1 + 2e) + ZyWk! ^ (^4) + 2e)'']ds^ 

_ ds .. ds“^ 

l^ ly 2e)'^dxdy 

+ [-Zyk- 0 . 2 s) 1/2 + 2 „k. # + 2 1^ ^ . 2 , |R k; f 
•'if "2wf ^^(1 .2.)]dxds 
+ [2xk- 0 . 2 s)’/ 2 - 2Wk; ^ . 2 f f - 2x f k; ^ 

^ 2y f k; ^ + 2W 1^ f ^ (1 * 2e)]dyds (81) 

Equations (80) and (81) define the metric coefficients g.. and 

■ J 

6.. for the undeformed and deformed beam, respectively, since the 

• J 

metric tensors g^j and are, by definition, given by 
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(d§)^ = g^.jdx’dx^ 
(dt')^ = G^jdx^dx^ 


Being that 


(ds')2 - (d§)^ = (G.j - g.j)dx’dx^ 



is a measure of the deformation, the covariant strain tensor is defined 
to be 


V. . — i(6* • ” g* •) (83) 

'ij 2' IJ 

However, the physical components of the infinitesimal strain tensor 
are needed for expressing the stress resultants in terms of the 
displacements. 

The normal strain components are defined by 



and the shearing strain components are given by 

where 0.. and 0'.. are the angles between the and x*^ 

I J I J 
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coordinate curves before and after deformation, respectively. These 
angles are related to the metric tensor components by 


cos 6, . = 


cos 61^ = 


i (86) 


J 


It is noted that in equations (84) and (86) and in the following 
expressions the repeated subscripts do not denote summation. Beginning 
with the definitions of equations (84) and (85) it is shown in reference 
8 and other texts that the normal and shearing strains are given by " 



Using equations (83) and (87) the above may be rewritten as 


®ii " 2 ^^ii ' ^ii^^^ii 


e . . “ [G 
ij ^ iJ 




(89) 


(90) 


As noted earlier, G. . and g. . are defined by equations (80) 

I J I J 

and (81). In the infinitesimal strain theory, products of deformation 
are neglected in arriving at equations (87) and (88). Thus, deformation 
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products will not be included in the expressions for the metric 
coefficients. . The linearized expressions for the curvature components 
kA, k', and k' are given by equations (62), (63), and (64). Also, 

A Jr . i 

e is defined by equation (21). These equations may be substituted 
into equation (81). to give (d§‘) In terms of the displacements a^, 
Oy, a^, and ({). Letting 

dx = dx^ 

2 

dy = dx 
dz = dx^ 


and comparing equations (80) and (81) with equations (82), the metric 
coefficients for the undeformed and deformed beam are found to be 


9n = 1 

922 = 1 

933 “ ’ ^ ''x - ''y) ‘‘x • •'z 



9i3 = -y •'z 


923 = ‘^z 
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^33 ” ^33 - xk^)(a^ + xk^^> + yk^cj) 

da da^ .2 ' 

^ ds V ds 

+ 2(-y + x^k^ + xyky)k^aj^ + 2(x + y^k^ + xyk^)k^Oy 


+ 2(x^ + y^ - xk^W - yk W)k 


y ' z 


®13 " ^13 - y(kx“x * ^“y ^ hI^ 

+ 0 + yk^ - xky) 1^ II + Wky ^ 

®23 ° ^23 * ’'('^x“x * '‘y“y * 

+ n + vk - xk ) — ^ - Wk ^ 
+ U + yK^ XK^J '^’^x ds 



From the above relations and equations (89) and (90) it is 
immediately seen that 



(93) 


It is assumed that the curvatures are small such that the product 
of a cross-sectional coordinate and a curvature component is much less 
than unity (i.e., x^ k. « 1). In other words, x and y are small 
compared to the radius of curvature and the angle obtained by 
multiplying the twist curvature k^ (radians per unit length of beam) 
by x (or y) is much less than one radian. With this assumption the 
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quantity appearing in equation (89) when i = 3 can be approxi- 

^33 

mated using the binomial theorem to give 


^ = 1 - 2(yk^ - xky) + 3(yk^ - xky)^ - (x^ + 

33 

Similarly, when e-j^ and e^^ a<^e determined, the inverse of the 
denominator of equation (90) is approximated by 


*■^11^33 ~ ^^13^ ^ ^ ~ ^^*^x ” ^^y^ ~ 2 ^ 


^^22^33 ” ^^23^ ^ - 1 - (ykj^ ~ ^^y^ "** ^^^x ^ ” 2 ^2 ^ 

(95) 

Besides giving more accurate strain expressions, the retaining of 

higher order terms allows the analysis, to be applied to problems having 

3 3 

larger curvatures. In the strain expressions terms of order x!j aj^ 

(or 0 ) will be kept. Equations (20) and equations (62) through (64) 

dot. 

indicate that k^. aj is of the same order of magnitude as . 

Thus, for terms involving derivatives of a , a , or <j», terms of order 

X y 

x? k? will be kept. Also, W will be taken to be 0(x?) since for all 

' I J I 

known solutions for the warping function, W is of the same order of 
magnitude as the product of the cross-sectional coordinates or smaller. 
Substituting equations (91), (92), (94), and (95) into equations (89) 
and (90) gives the strains ®2Z^^®33^’ ®xz^"®13^’ ®yz^“®23^ 
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= [-y - xyky + 2y^kj^ + - x^yk^ - 3y^k^ 

\ 

+ 2x^kj^k - 2x^yk^ + 4xy^k^k^ + y(x^ + 

+ [x - xykjj + 2x^k + y^k + xy^k^ + 3x^k^ 

- 2y^kj^k + 2xy^ky - 4x^ykj^ky - x(x^ + y^)k^]kj,Oy 

+ [1 + xky - ykj^ + x^k^ - 2xyk^ky + y\\ - (x^ + y^)k^ 

+ x\® - 3x^yk^k^ + 3xy^k^ky - y^k^ + 3(x^ + y^)(yk^ - xky)k^]a^ 
+ [xkj^ + yky + x^kjjky + xyky - xyk^ - ■*■ ’‘^'^x'^y * ’’^^‘‘y 

■ i • * 

- 2x^yk^k - 2xy^kj^ky + xy^k^ * /‘‘x'^y ' ('‘‘'x 

- x[l - (yk^ - xky) + (yk^ - xk^)^ - (x^ + y^)k^] 

- y[l - (yl^x * ’‘‘‘y^ * ^^*^x ' * y 3i^ 

■■ . • . _ . • ■ 

' j 2 ' 

+ [{1 - 2yk^ + 2xky)(x2 + /) - (xk^ + yk^)W]k^ ^ + W (96) 
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®xz = -ytl'x - (y"x - '‘^)''x yl'z + (yl'x - '“'y)^''x 

+ (- I x^kj^ + 2xyky - 3y^kj^)k^]aj^ 

- y[ky - (ykx - xky)ky - xk^ + (yk^^ - xk^j^k^ 

+ (- I x^ky + 2xyk^ - y^ky)k^]ay 
+ yk^n - 2(ykj^ - xky) + 3(ykj^ ' ‘ (f + y^)k^]a^, 

+ yk^Cxkj^ + yky + 2{x^ - y^)kj^ky + 2xy(k^ - k^)](f. 

- xyk^Cl - 2(yk^ - xky)] ^ - y^k^p - 2{yk^ - xk^)] ^ 

- y[l - yk^ + xky + (ykjj - xky)^ - (| x^ + y^)k|] ^ 

" 0 - i f ^ + (1 - yk^ + xky)Wky ^ + yWk^ ^ (97) 
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®yz = 

+ (- |- + 2xyky - x\j^)k^]a^ 

+ x[ky - (yk^ - xky)ky - xk^ + (yk^^ - xk^)^ky 

+ (- f y\ + 2xyk^ - 3x2ky)k^]ay 

- xk^[l - 2(ykj^ - xky) + 3(ykj^ - xk^)^ - (x^ + | y^)k^]a^ 

- xk^[xkjj + yky + 2(x^ - y^)k^k^ + 2xy(k^ - k^)]4> 

9 <ict„ da, 

+ X%[1 - 2(yk^ - xky)] + xyk^El - 2(yk,^ - xk^)] g^L 

+ x[l - ykjj + xk^ + (ykj^ - xk^)^ - (x^ + | ^ 

+ 0 - j y^l«z)|y ^ - n - yk^ + xky)Wk^ ^ - xWk^ ^ (98) 

In the above strain-displacement relations terms have been collected 
on each variable. If the curvatures are very small it is obvious that 
the majority of the terms are negligible. When retaining only the 
largest terms associated with each variable the strain expressions 
are approximated by 
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da 

«2z “ “z - >'(-''z“y - Vz - ‘'x* di-) 


+ y(-k a - ka + k^“ (x^ ■*■ y^)k_ ^ W 

jv ■'^“x X z y^ ds ' ' ' z ds 


( 99 ) 


®xz = Vy ■ '"z“z ^ ^ lx ^ 


da « da j2. 

- "‘■''“z ar ■ y "z di^ * y'^z^’^'^x ■' y"z ” ^ 


( 100 ) 


®yz ° "'(''x“x •" ^“y - ‘'z“z - I 7 ^ 


2.. <*“x 


da 


d^<t) 


^ dT ^ ^^*^2 dT - ^ ^ 


(ion 


D. Discussion of Straln-Displacement Relations 
In the expression for e^^ as given by equation (99), the first 
term is the linearized form of the strain of the elastic axis that was 
given in the development of the curvatures. The quantities inside the 
first two sets of parentheses are equal to (k^ - k^) and (k^ - k^), 
respectively, and these are the changes in the curvatures of the elas- 
tic axis due to deformation. The term (x + y )k^ d4)/ds is due to 
the initial twist of the beam. This term has been obtained previously 
by Houbolt and Brooks (see reference 14) in their analysis of a 
straight pretwisted propeller. The last term in equation (99) arises 
because uniform torsion has not been assumed. 
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The quantity enclosed in the first set of parentheses in each of 

equation (100) and (101) is equal to the change in the torsional 

3W d(b j 3W d^i 

curvature component (i.e«, - k^). The terms ^ and ^ 

contribute to the Saint-Venant torsional stiffness constant that appears 

in the torsional stress resultant. Based on the previous observation 

da^ 


that k. a- may be of the same order of magnitude as and the 

assumption that k. x. « 1, the remaining four terms in each of 

1 J 

equations (100) and (101) may be smaller than those mentioned above. 


However, these terms cannot be assumed negligible since there may be 

instances where, for example, may be zero while its derivative is 

nonzero. The last four terms in each of equations (100) and (101) i 

appear because of the beam's initial - twist curvature. 

It is noted that since shear deformation was neglected the strains 

p and e do not reflect the strain arising from transverse shear, 
'^xz '"yz 

The expressions for e^^ and are used to determine the torsional 

moment stress resultant. The omission of the strain due to transverse 
shear does not affect this stress resultant since moments are taken 

about the shear center of the cross-section. 

In the development of the strain relations it was assumed that the 
warping displacement was proportional to the Saint-Venant warping 
function W. The solution for the Saint-Venant warping function of a 
cross-section assumes that the torsion is uniform (constant torque, 
constant d<{)/ds). Since this analysis is not restricted to uniform 
torsion there is an inconsistency. However, for many problems in which 
the torsion is nonuniform, it is common practice to determine the 
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torsional stiffness using the Saint-Venant theory. This is especially 

true for natural vibration problems. The effect of nonuniform torsion 

on the torsional stress resultant becomes important for flanged thin- 

wall open sections such as I-beams. For such beams the theory has been 

developed to correct the stress-resul tant expression for nonuniform 

torsion (see reference 15). For beams which are not thin-walled, a 

general theory for nonuniform torsion has not yet been developed. For 

this analysis, it will be assumed that the torques acting on the beam 

are distributed loadings (such as inertial loads) and that the derivar 

tives of d<J)/ds are small such that the effects of nonuniform torsion 

2 2 

are negligible. With this assumption, the terms proportional to d <|i/ds 
in the strain expressions may be dropped. However, for the examples 
presented herein these terms have no effect on the stress resultants 
since the cross-section is assumed to be doubly symmetric. 

The strain-displacement relations have been developed for the pur- 
pose of obtaining expressions for the stress resultants in terms of 
the displacements. In the classical theory for curved space beams the 
strain-displacement relations are not derived. The stress-resultants 
are assumed to be equal to the changes in the curvatures (without the 
extensional effects included) times a cross-sectional stiffness 
property. The importance of using the strain-displacement relations 
instead of the classical approach will be seen when the stress resul- 
tants are compared since the classical approach can lead to considerable 


error. 
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It was noted earlier that the accuracy of the strain expressions 
depends upon the number of terms kept 1n the binomial expansion of the 
denominators of ec^uations (89) and (90), The strain e^^ as given by 
equation (96) can be compared with the Winkler theory for curved beams 
The Winkler theory, which is valid for large curvatures, .applies to 


untwisted beams lying in a plane, having a constant radius of curvature, 
and with the deformations being in-plane. Thus, considering a beam of 
constant curvature which lies in and is displaced in the x, z plane 
such that ky = 1/R and = k^ = Y = v ^ equation (96) 

gives 


zz 


= 0 + 


2L 

R 




- x(l +1 


2 da 
X \ 

^ ds 


'I 


(102) 


Substituting equations (20) for and yields 


zz 


- dw /•, . X + ^ + 2L-) y. - x(l + w 


3' R 


x^t d^u 

^ 7^ ds2 


(103) 


It is noted that according to the sign convention established by 
figure 2 and equations (1), x and u are positive in the direction 

toward the center of curvature when y = 0. 

In reference 16 Langhaar gives the following expression (notation 
and signs have altered to conform to that presented herein) for the 
strain according to the Winkler theory 


2 2 

dw u/R ^ X d u/ds 
Sz “ ds " 1 - x/R “ 1 - x/R 


( 104 ) 
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The above gives a parabolic distribution of strain over the cross- 
section instead of the familiar linear distribution associated with the 
bending of straight beams. This distribution is brought about because 
equation (104) was obtained by taking into account the differences in 
the initial lengths of the longitudinal fibers in a differential segment 
of the beam. The strain distribution given in equation (103) appears 
because approximations were made for the complicated denominators of 

the strain expressions for space beams. 

From the binomial theorem it is seen that equations (103) and (104) 
are essentially the same when x is small compared to R. In fact, 
since the higher order terms were kept, the strain expressions can be 
applied to beams having moderately large curvatures. For example, 
suppose that x = 0.2 R. Then, equation (103) gives 

1-248^- 0.248 

ds 

From equation (104), which Is the expression for large curvatures, it 
is found for x/R = 0.2 that 

®zz ^ dF ~ 1«250 ^ - 0.250 R 

ds 

If X = 0.2 R is largest cross-sectional coordinate of a point in the 
beam, then it is seen by comparing the two above results that the 
largest error in the strain as given by equation (103) is less than 
one percent (the same error is obtained at x = -0.2 R). Equation (103) 
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is extremely accurate for x/R £0.1 and does not give grossly 
inaccurate results until x Is about one-half of R. , 

Equation (99) gave the expression for e^^ with higher order terms 
dropped. For the example being considered it reduces to 


®Z 2 ds ^ R ^ 


At X = 0. 2 R the above becomes 


= ^ - 1.2 §• - 0.2 ^ 
ds R jj2 


When compared to the result of the Winkler theory, the last term is 
seen to be 20% in error. 

In the development of the stress resultants the more>accurate 
strain expressions containing the higher order terms will be used (i.e., 
equations (96), (97), and (98)). 

It is recalled that in the classical theory the elastic axis is 
assumed to be inextensible such that = 0. From equation (102) it 
is observed that the extensional strain ct^ of the elastic axis must 
be included. Otherwise, the strain-displacement relation bears no 
resemblance to what it should be. 

In reference 17 Oden gives an expression for the strain of a 

curved planar beam of radius R. His result is different from that 

1 ' 

given by Langhaar and from that presented herein. Oden arrives at the 

I 

total strain by summing strains that arise from threei types of 
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displacement (I.e., extensional, radial, and rotational)- The strain 
due to the radial displacement is given as 

\ _ -y V 
" 1 ~ y/R j^2 

where v is the inward radial displacement- The strain due to a radial 
displacement is actually the second term of equation (104). By 
considering the strain due to the radial displacement of a circular 
ring it is obvious that the above result is in error. The error of 
reference 17 was caused by taking the final, length of an incremental 
fiber at a distance y from the centroid to be As[l - y(R - v)] 
instead of As[l - (y + v)/R] as it should be. 

E. Development of Stress Resultants 
The stress resultants are determined in the usual manner by 
summing over the cross-section the stresses that are acting on it. But 
first it is noted that in general the x, y, and z curvilinear 
coordinates (s and z are used interchangeably) and not orthogonal. 
When the beam is twisted (k^, f 0) the coordinate lines are orthogonal 
only at points along the elastic axis. The x and y coordinate lines 
are orthogonal at every point in the cross-section before deformation 
(g-j 2 = 0) and after deformation (G ^2 ^ ^ when neglecting product of 
deformation). But the z coordinate line passing through a point in 

I 

the cross-section off the elastic axis is not normal to the xy plane 
(9^3 ^ 0» 923 ^ 0» ^ ^23 ^ can be visualized for a 



65 


twisted beam since the z coordinate line is defined by the inter- 
section of the coordinate surfaces x = constant and y = constant. 

Figure 3 shows the stress acting at point P in the cross- 

section and in the direction of the z coordinate line passing through 
that point. Also shown are the components of in the directions 

of e‘, e', and e*. The components of o_ in the e', e‘, and e* 

-wZ ^jr -srZ ZZ -vX 

directions contribute to the torsional moment stress resultant. This 


effect, which occurs in twisted beams, is discussed by Goodier in 
reference 15. The bending moment stress resultants about the x and 
y axes and the axial force stress resultant are determined from the 


component of in the e^ direction. 

The component of a in the e* direction, for example, is a _ • 

ZZ ZZ 

cos 0 -J 2 where cos 0^^ is defined in terms of the metric coefficients 


G.. by the second of equations (86), However, since products of 
deformations will not be Included in the stress resultant expressions, 


the reference of the undeformed beam will be used. 


If the subscripts 


Xq, y^, and are used to denote the components of at point P 

in the e^, e , and e directions, respectively, then 


^° 2 z\ = °zz ®13 = “zz^ia/*^ 

“ °ZZ ®23 ' °zz^Z3^'^ 

t 

since g^^ = 9^2 “ 1- Letting 1^ be the direction cosine of the angle 
between o and e , the orthogonality relation . 




s 



Figure 3.- The stress at a point in the cross-section. 
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cos^ + cos^ 6^3 + 1^ = 1 


gives 


? 2 1/2 

(y 

' 933 ' 


Thus. 


= ®zz ” 
‘o ■ 


('- 


>,2 + „2 
^13 ^ ^23 

% 


1/2 


(106) 


From these components of and the shear stresses acting in 

the X, y plane, the stress resultants are given by the following 
integrals. 

o o 1/2 






913 * 923 \ 
933 / 

, 2^2, 

ez^dA 

1/2 

(107) 

9-|3 + 923 

^ 933 i 

) ^zz^*'' 

(108) 

^9i3 , *923 \ 
^ ^3' 

ez,dA 

(109) 

J- ^ .2 . 1/2 


yi3 ^ ^23 
933 / 

«zz 

(110) 


Using the binomial theorem, it is found from equations (91) and 


(96) that 
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(l ^13 
' ■ 933 / 


zz 


= [-y - xyky + 2y^kj^ + x^k^ - x^yk^ - 3y^k^- 

+ 2x^kjjky - 2xVx + ^xy^k^ky + | y(x^ + y^)k^]kj,a^ 

+ [x - xykjj + 2x^ky + y^ky + xy^k^ + 

- 2y^kyky + 2xy^ky - 4xV^ky - | x(x^ + y^)k^]k^a^ 
+ [1 + xky - yky + x^ky - 2xykj^ky + y^k^ 

- I (x^ + y^)k^ + x^k^ - 3xV^ky + ^xy^k^ky 

- y^k^ + I (x^ + y^)(yk^ - xky)k2]a- 

+ [xkjj + yky + x^kjjky + xyky - xyk^ - y^k^^k^ 

+ x^kjjk^ + x^yk^ - 2x^yk^ky - 2xy^kyk^ 

+ xy^k^ + y^k^k - | (x^ + y^)(xkj^ + yk )k^]((> 


- x[l - ykj^ + xk + (yk^ ' xk )^ - | (x^ + y^)k^] 


da 


- y[l - yk^ + xky + (yk^ - xky)^ - | (x^ + y^)k^] ^ 

2 

+ C(x^ + y^)(l - 2yk^ + 2xky) - (xk^ + yky)W]kj, ^ + W ^ 


Oil) 
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^ ” . 

= {-y - 2xyky + 3y^kj^ + x^kj^)(x^ + y^lk^a^ 

+ (x - 2xykjj + 3x^ky + y^ky)(x^ + 

+ [1 + 2xky - 2ykj^ + 3x^k^ - Sxykj^k^ 

+ 3y^k^ - I (x^ + y^)k^](x^ + y^)k^a^ 

■ 1 

+ (xkj^ + ykyHI - 2ykjj + 2xky)(x^ + /)k^4: 

- x(l - Zy\ + 2xky)(x^ + y^)k^ ^ 

O P 

- y(l - 2ykj^ + 2xky)(x^ + y )k^ ^ 

+ (x^ + y^)^ k^ # + (x^ + y^)Wk ^ (112) 

z as z jjjZ 

Froni the above equations it Is seen that the stress resultant 
expressions will be rather lengthy when the above and the shear strains 
are substituted. Solely for the purpose of reducing the size of the 
stress resultant expressions it is assumed that the cross-section has 
two axes of symmetry. Thus, the following Integrals are zero. 


Ill I 
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f xydA = f x^dA = f 
JA JlK Jp, 


f xdA = f ydA = f xydA = f x^dA = f y^dA = 0 
M Ja 




r x^ydA = f xy^dA = f x^ydA = f xy'^dA = 0 

Ja Ja Ja Ja 


) (113) 


J 


Also, the warping function W(x, y) is an odd function of both x 

' . CJ 

and y since the warping displacement must be antisymmetric with 
respect to both the x and y axes for doubly symmetric cross- 
sections, Thus, o 


f WdA = f xWdA = f yWdA ■= f x^WdA = f 
•'A /A Ja Ja y 


y‘WdA = 0 


(114) 


Even if the cross-section Is not symmetric the second and third 
integrals above are zero since x and y are measured from the shear 
center (or center of twist). This was proven by Goodier in reference 
18, 

It is assumed that the elastic modulus E and the shear modulus 
6 are constant over the cross-section. The integrals which appear in 
the stress resultants are denoted as follows 



1 
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A = f 

JP 

dA 


>1 

'xx = 

f y^dA 

-'A 




f x^dA 

•'A 



'p = J 

(x^ + y^)dA 
A 




(x2 + y2 + X f 




x^dA 
A . 




I* y^^dA 

A 


> 

B3 = J 

x^y^dA 

A 



«4=J 

^ x^(x^ + y^)dA 
A 



^5-J 

y2(x2 + y^)dA 

A 




(x2 ^ y2j2j^ 

A 



«7 = J 

xyWdA 

A 



CO 

II 

i " 



®9 = J 



j 
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In the above, J is the well-known Saint-Venant torsional stiffness , 
constant. 

Substituting equations (97), (98), (111), and (112) into equations 
(107) through (110) and applying equations (113) through (115) and 
replacing a_ by its definition 


“z = ^ - 


yields the following expressions for the stress resultants. 


”x = * 2B3k2 - I B5k2)k^a^ 


2(B2 + 2B3)k^kyk^cXy 


(^xx" V5"¥M V>x^ 




- ('xx - I 


* (‘xx " ^3^5 - + B2k2 - | B3k2)ky$ 

9 2 3 2 

^ 2B3kx^ dT - ('xx * * h< - f d/ 


- (2^5 ^ B,)k^k^ 


(116) 
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My = E[-2(B^ + 2B3)kj^y^a^ 

- < ^ ^'3'^y - I «4^>z“y 

- (Iyy-^BikJ + 3B3k2-fB,k>yf 
Ml,y + B^kJ*3B3k2-|B,k2)k2u 

^^yy * ^I'^y * ^®3*‘x ‘ I Vz^*'x''y'' 

-{Iyy-.B,k5-2B3k2 + B3k2-|B4kX^ 

9 9 9 

" (Iyy-^'lV®3'<M 

riA 

-^Wy^- (2B4 - B7>kyk^ ^] (lip 

«z = + 3B2 f 5B3)k2 . I OB5 + B,)k2]k^a^ 

+ [Ip + B^k^ + Bgk^ - (3B^ + B3 + 5B3)k^ + I {3B^ + Bglk^lkyay 

+ [-Ip - BB^kJ - BBgk^ + (B3 + B^kf 

-|(Ip*3B,k2.3B5k2-|Bgk2)]k^^ 

- [-Ip - 3B^kJ - 3Bgk2 + (B 3 + Bglk^ 

* I + 3Bgk^ - I Bgk^)]kyk^u 

+ [-Ip - 3B^ky - 3Bjkj^ + (Bj + Bg)k^ 

+ |(Ip^^V5^3B5k2.|Bgk2)]k^k^v 

+ 2 [B 2 - B^ + I (B 4 - Bgnk^kyk^* 
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da 


+ 2{1 - - 2(1 - ^5 


da 


+ [J + (B^ - B^jk^ + (Bg + Byjk^ - (Bg + Sg - ^ 8g + |- Sg)k2 


E B .2-, # 
^ G ®6'^zJ ds 




(118) 


>^z " Ixx^'^z“x " (Ip " VVz“y 

+ fA + I + I - — I ^ 

^ yy y xx\ i ds 
- " Vy " 'xx>^x -fVzV 

^ " ^x*^x - I V\K^ 


da 


da 


1. + T k 

yy ' ^x'"xV ■ Vy"y ds' ‘ *xx^ ds ^ ds-* 


+ (I.,., - Ivv)k„k * - I.,„k., 1 


(119) 


F. Discussion of Stress Resultants 
Equations (116) through (119) give the stress resultants in terms 

of the displacements. These expressions include extensional effects 

■ *1 

which have not been taken into account previously in stress resultants 
for space beams. Also, the expressions contain terms whiph allow the 
equations to.be applied to beams having curvatures larger, than . 
would otherwise be permitted. 

It was noted earlier that in the classical approach the stress re- 
sultants are not derived from strain-displacement relations but are as- 


sumed to be proportional to the changes in the curvatures. Thus, accord 
ing to the classical theory, which also assumes the elastic axis to be 
inextensible (a^ = 0), the three moment stress resultants are given by 
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da 

M = El „(-k a + k 4 - 3-^) 
X XX ^ z X ds ' 


da 

M = El (-k a - k d) + 

y yy z y ^ ds ' 


) ( 120 ) 


M = 6J(k a + k a + ^) 
z ' X X y y ds/ 


J 


It is obvious that the more rigorous approach for determining the 
stress resultant yields many additional terms. In equations (116) 
through (119) the dw/ds, u, and v terms appear because of extension 
of the elastic axis. The terms Involving the cross-sectional 
constants arise because of the more accurate representation of the strain 
distribution over the cross-section. The first two of equations (120) 
essentially assume a linear strain distribution. When the cross- 
sectional dimensions are very, very small compared to the radius of 
curvature (i.e., x^ kj is much, much less than unity) most of the terms 
In the stress resultants which contain a constant can be neglected. 

A comparison of the third of equations (120) with equation (118) 
points out that a very common mistake, which has not been noted before, 
has been made in previous analyses that dealt with torsional and 
out-of-plane deformations. The comparison shows that in the third of 
equations (120) only the displacement d4>/ds should be multiplied by 
the torsional stiffness 0. The variables a and a appearing in 

^ y 

this equation should be multiplied by the polar moment of inertia I . 

P 

Unless the cross-section is circular, the use of the stress resultant 
Mj, as given by the classical approach can lead to considerable error 
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since Ip may be. many times larger than J when the cross-section is 
not compact. Consider, for example, a rectangular cross-section. The 
ratio of I to J is given below for various ratios of the cross- 

K 

sectional dimensions a and b. 


a/b 

ip/j 

1 

1.19 

2 

1.82 

3 

3.17 

4 

5.04 

5 

7.44 

10 

26.98 


The easily made mistake described above has often been made in ^ \ 
papers pertaining to torsional and out-of -plane vibrations and out-of- 
plane buckling (for example, see references 12, 19, 20, 21, 22, and 
23). In the classical formulation the fact that the a and a 

^ y 

terms should be multiplied by Ip instead of J is not at all pbvious 
This error shows the need for determining the stress resultants from 
strain-displacement relations. 

In reference 15 Goodier showed that for a straight twisted beam 
the pretwist contributes to the torsional stiffness since the displace- 
ment d(j)/ds causes a longitudinal stress which opposes the motion. In 
the notation of this paper, Goodier gives the stress resultant as 


<= (GO + EB,k^) ^ + Ei„k ^ 
z ' o z' ds p z ds 
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These terms, which appear in equation (118), are shown to be significant 
for slender cross-section, such as propellers, even when times the 
largest cross-sectional coordinate is as small as 0.1. In the stress 
resultant development, the inclination of causes terms other than 

those above to appear in since the treated beam also has the 
curvature components k and k . 

X y 

In reference 24 Den Hartog gives an interesting observation that 
he made in his experiments with bending of pretwisted straight beams 
and challenges the reader for an explanation. It was found that the 
bending stiffness of doubly symmetric beams having equal cross-sectional 
moments of inertia is reduced when the beam is pretwisted. The cruci- 
form cross-section is cited as an example. According to the conven- 
tional theory for straight pretwisted beams (for example, see 
reference 25), the pretwist has no effect on the bending stiffness (nor 
on the bending displacements) when the moments of. inertia are equal. 

In reference 26 the reduction in stiffness is attributed to distortion 
of the cross-section and anti clastic effects. Equations (116) and (117) 
show that there is also a reduction in the bending stiffness due to the 
pretwist itself. For a straight pretwisted beam having only the dis- 
placements u and V, the equations give 

»x = -Edxx - f Vz) 

MyVE(Iyy-|B,k2)^ 
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It can be shown that for the cruciform cross-section, and Be are 
equal, as well as and and do not depend upon the orientation 

chosen for the x and y axes. Thus, it is seen that the pretwist 
does indeed affect the bending stiffness even when the moments of inertia 
are equal. 



The equilibrium equations are derived by considering the forces 


and moments acting on an incremental segment of the beam. The equations 

j j 

are obtained using a very expedient vector approach. 


The internal moments (or moment stress resultants) acting in the 
?x* Sy’ ?z ^■''^^ctions were previously denoted by M^, M , and 

•J 

M^, respectively. Similarly, let F^, and F^ denote the internal 


forces acting in these directions. Also, let p , p , p , q , q , and 

x y 2 . X y 

q^ be the applied distributed forces and moments acting in the e', e’. 


~x’ ^y 


and e^ directions per unit length of the extended rod. From these 


components mentioned above the following vectors may be defined. 


* Vy + ”z?z. 

P = Px!i * PyP; ^ Pz!z 

3 = Px!; * ^ 
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Figure 4 shows an incremental segment of the beam of length As'. 
Thus, the applied force and applied moment vectors acting on the incre-- 
ment are p As' and q As'. The changes in the internal force and 
internal moment vectors from one side of increment to other are denoted 
by AF and AM as shown. To establish the sign convention, consider 
an imaginary cross-sectional cut of the beam taken at s' = sj. The 
sign convention is chosen such that the internal force and moment 
components acting on the cross-sectional face of the material on the 
s' < s^ side of the cut are positive in the e^, e^, and e^ 
directions. 

The summation of forces and the summation of moments about the | 

i. ^ 

left hand end of the Increment give 

F + AF + p As ' - F = 0 
M + AM + q As’ + (As'ep x (F + AF) 

+ (|- As'ep X (p As') - M = 0 

(f 

Dividing each of the above equations by As' and taking the limit as 
As' approaches zero gives the equilibrium equations in vector form as 

dF 

g|r+P = 0 (122) 

dM 

e; X F + q = 0 


(123) 
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Substituting equations (121) into the above, performing the indicated 
operations, and substituting equations (32) for the derivatives of the 
unit vectors yields 


& - * Px = 0 

:t-+ + k'F - k'F + p » 0 
ds^ z X X z 


dF. 


- k.’F. + + P, = 0 


^ Vx Vy 


> (124) 


dM. 


- klM,, + k‘M, - F + q = 0 


ds‘ ’'z’V Vz y 


dM 

^ + k;M^ . = 0 




> (125) 


The above equilibrium equations are the same as those given by 
Love in reference 1 except that ds‘ appears in the above derivatives 


instead of ds since extensional behavior is considered. Replacing 
ds* by 


ds' = (1 + 2 ds 


in the above gives equations of the type 
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dF 

+ 0 + + k'F^ + Px) = 0 (’26) 

As noted before, is the force acting in the e^ direction per 
unit length of the extended rod. Let dP^ be the force acting on ds'. 
That is, 

Px “ d?*” 

Therefore, 

0 . p, ■ (1 . 2.)'« M 

' . ■' i 

. dP 

■ ' s X 

' ' ■ ■ ” ds 

Let dP /ds, which is the force acting in the e* direction per unit 

X , 

length of the undeformed beam,' be denoted by p^. Thus, 

(1 + 2e)’/2 (127) 

There are similar relations for p^, p^» g^i and q^. Hence, the 
equilibrium equations as typified by equation (126) take the form 

+ (1 + 2e)’/2(-k-Fy + k'F^) + = 0 


(128) 
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Suppose that is the inertial load given by (-m*u) where m* 
is the distributed mass of the extended beam. Then, 

(1 + = -(1 + m‘u 

Let dM be the mass of the beam in the differential lengths ds and 
ds*. Using equation (127), the above becomes 

p^ = .(, .2e)V2^a 


= -mu 

where m is the distributed mass of the undeformed beam. . Therefore, 
when considering inertia loads the distributed mass in the p , p , 

^ y 

etc., expressions is that of the undeformed beam. 

In this analysis only linear solutions are being considered. For 
free vibration problems the internal forces and moments are proportional 
to the displacement variables. In the and k' expressions, 

^ jr Z 

the initial curvatures are the only constant terms. Thus, from equation 
(128) and those similar to it, the linear equilibrium equations for 
natural vibrations are found to be 
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dF. 


— - - k F 
ds y 


— jL + r 

ds 


dF. 


- k..F. 


ds Vx 


+ k F + p = 0 
y z 

- k F + p - 0 

X 2 

+ k F + p = 0 
X y ^2 


> ( 129 ) 


ds 


. kyM^ 




= 0 




di^ ^ ‘‘z”x 


dM^ 

ds 


kyMx 


- Vz " 'x ; s = ° 

+ k^Hy + = 0 


> (130) 


H. Surrmary of Governing Differential 
Equations for Natural Vibrations 

All of the equations needed for solution have now been developed. 

It is planned for the natural vibration solution to be obtained using 

a transfer matrix technique. Thus, the first-order equations that have 

been derived will not be combined to form higher-order equations. The 

problem has twelve variables; u, v, w, (j), oiy, F^, F^, 

and F^. The twelve equations are: the first two of equations (20) 

which relate a and a to u, v. and w; equations (116) through 
X y 

(119) which are the four stress resultant relations; and equations 
(129) and (130) which are the six equilibrium equations. 

The first two of equations (20) may be rewritten as 
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^ = k^v - kyW + 


^ = -k^u + k^w + ay 


) (131) 


The four stress resultant relations can be written in matrix form 


as 


M 


r ^ 
u 


ds 


> = ta] / ^ I + [b] ) 


2 

^ J 


a. 


a. 


y 

V J 


di 

ds 


dot^ 

dT 


da 

y 

ds 

J 


(132) 


where the matrix elements a^j and 
through (119). It can be shown that 
singular. Thus, the above becomes 


b^j are defined by equations (116) 
the matrix [b] is always non- 
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r N 
dw 
ds 

' 

u 

' 


ds 

> = - [b]-’[a] ( 

V 

> + [b]'’ { 


da 
X 

ds 

“x 

, 





da 

di^ 

J 

■ ' 

“y 
V. J 




The distributed applied loads p , q , etc,, appearing In the 
equilibrium equations consist of inertial loads and elastic foundation' 
forces and moments. Rotary Inertia and elastic foundation effects are 
included in the equations since they do not complicate the transfer 
matrix solution In any way. 

The distributed mass m is equal to pA where p is the mass 
density and A Is the cross-sectional area of the beam. The reversed 
effective inertial force In the x-direction, for example, Is (-pA 3^ u/3t^), 
Thus, for harmonic motion 


2 I, 

= u) pAu - c^u 
2 

Py = U) pAV - CyV 
2 

P^ = 0} pAw - c^w 


>( 134 ) 


J 
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where w is the natural vibration frequency and the c^. are the elastic 
foundation constants for translational displacements. 

In the development of the curvature relations the Eulerian rota- 
tions about the x, y, and z axes were denoted by 0, and 4>, 
respectively. The rotations 0 and are given in terms of the elastic 
axis displacements by equations (45). In the linearized analysis e 
and are small and are approximated by 



' “x 


The rotary inertias about the x, y, and z axes, respectively, are 

1 ■ »‘yy , 

¥ 

f p(x^ + y^)dA = pi 

A P 

D 

Also, it is assumed that the elastic foundation can resist rotations 
about the x, y, and z axes, and these foundation constants are 
denoted by d^. Thus, using the above relations for the rotations Q 
and Ip, the applied moments are .given by 
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Q - -w dI a . + d a 
^ XX y X y 


Q = w^dI a - d ex 

Hy " ^ yy X y x 


2 

- w plp<|) - d^cj) 




) ( 135 ) 


Substituting equations (134) and (135) into the equilibrium rela- 
tions as given by equations (129) and (130) and rearranging yields 


J = XA - ^ - ‘‘x>“y 


dM 9 

^ =-kM +kM -F - (w^pl - d )a 
ds z X X z X ^ ^ yy y^ x 


dM 

3T ' •'yMx - ^«y ' ' “zJ* 


2 

dT ' ''z'^y * '"y'^z ' P* ' Px>“ 


‘‘''v 2 

= -kF +kF - (w^pA - c )v 
ds 2 x x 2 ^ ^ y' 


&Wx- 


) (136) 






) (137) 


Equations (131), (133), (136), and (137) are the twelve first- 
order differential equations which describe the natural vibration 
motion of space beams. The motion consists of bending in two 
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directions, extension, and torsion. The degree to which these motions 
are coupled depends upon the application. It Is noted that the usual 
beam bending-torsion coupling arising from the noncoincidence of the 
elastic and centroldal axes does not appear because of the assumption of 
double symmetry for the cross-section. 

I. Numerical Solution Method 

As noted earlier the solutions for the natural vibration character- 
istics win be obtained using a transfer matrix method. This method, 
which Is based on second -order Runge-Kutta Integration, does not 
discretize (I.e., lump masses and assume constant elastic properties 
for a beam segment) the beam as most other transfer matrix methods do. 
Instead, distributed mass values and elastic properties at selected 
stations are used. The method assumes that the mass and all other 
properties vary linearly from station to station. Also, any discon- 
tinuities In the beam's properties are conveniently handled by assign- 
ing two stations to the point at which the discontinuity occurs. In 
addition, the method can be applied to rings as well as to beams. 

Matrix differential eouatlon .- Let {Y} be the column vector 
consisting of the twelve variables. That Is, 
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{Y(s)} = 




(138) 


Also, let {Y'} be a similar vector defined by the derivatives of the 

i 

variables. The set of differential equations given by equations (131), 
(133), (136), and (137) may be expressed in matrix notation as 


{Y*(s)} = [H(s, (xi)]{Y(s)} 


(139) 


Subscripts will be used to denote stations. Thus, at station i 
the above is written as 


{Vp = [Hp{Yp 


(140) 


Development of transfer matrix .- The value of {Y} at station 
(i + 1) may be approximated as follows. 


{Y 


i^l} = {Yi}+^{tYp-MYJ^l}} (141) 


I 
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Equation (139) may be written at station (i + 1) and substituted, along 
with equation (140), into equation (141) to give 

Now, the on the right hand side of the above equation may be 

approximated as 

(Vi) = {Y^} + ASi{vp 

= {Y^) + As^tH^]{Y^} , 

Substituting into equation (142) gives 

= {Y,} + ^ |[H^]{Y^} + [H^+i]{(Y^} + As^[H^]{Y.)}| (U3) 

The above may be rewritten as 

= [Ai]{Y^> (144) 


where 


[A^] = [I] + ^ As^J[H^] + [^5+]]] r ^^i 
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The above matrix, [A^], is the desired transfer matrix* It is used to 


relate the variables defined by {Y} evaluated at station (i 
those at station i. 

+ T) to 

Beam solution.- Lettinq i = 0 and i = n denote the stations at 
the ends of the beam, {Y^> may be related to {Y^} by repeated use of 
equation (144). Hence, 

{Y„} ' [B]{Y„} 

046) 

where 


[B] = . . . [AglLA^lCA^] 

(147) 

The six homogeneous boundary conditions at each end of 
may be expressed in matrix form as 

the beam 

[C]{Y^} = {0} 

(148) 

[D]{Y(j) » {0} . 

(149) 


12 matrices. Equation (146) may be 
(148) and the result combined with equation 

[U]{Yq> = .{0} (150) 


where [C] and [D] are 6 x 
substituted into equation 
(149) to give 
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where 


[U] = 


[C][B] 


[D] 


(151) 


For a non-trivial solution to equation (150) the determinant of 
[U] must be zero. In the usual fashion for the transfer matrix method, 
the values of ui which make 


|[U]|=0 (152) 

i 

are found by iteration. For each natural frequency the solution for 
{Y..} at every station is obtained by first setting one element of {Y^} 
equal to unity and solving equation (150) for the remaining elements of 
{Y^}. Then the {Y..} are given by equation (144). Also, the derivatives, 
{Yj}, at each station may be computed from equation (140). 

Ring Solution .- For rings and other curved beams which do not have 
endpoints the "boundary conditions" are applied in a different manner 
than that above. One arbitrary point on the ring is chosen to be both 
station i » 0 and i = n. The variables of {Y} at these stations 
are related by 

{Y„} = [E]{Yjj} (153) 
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where E is a 12 x 12 matrix. When the total twist of a ring is zero 
or a multiple of 2 it radians the matrix [E] is the identity matrix. If 
the total twist of a ring is an odd numbered multiple of tr radians 
then the matrix [E] is the negative of the identity matrix. Substitut- 
ing equation (146) into the above gives 

[B]{Yq}.= [E]{Y^}. 

[[B] - [E]j{Yp} = {0} (154) 

other than the trivial solution, 

|[B] - [E]!' = 0 ; . . (155) 

As before, the values of m which make the above determinant equal to 
zero are found by iteration. The modal characteristics are determined 
using equations (154), (144), and (140). 


Or, 


For solutions 



VIII. NUMERICAL EXAMPLES 


Two numerical examples are given which utilize the deformation 
theory and numerical method presented in this dissertation, A twisted 
curved beam and a twisted ring are considered. It is thought that 
natural vibration solutions have not been obtained previously for beams 
which are both curved and twisted. 

Twisted curved beam.> The natural vibration characteristics were 
computed for a beam whose elastic axis forms a semi-circle as shown in 
figure 5. Both ends of the beam have cantilever boundary conditions 
(u-v-w-(|i = a^ = ay = 0). The radius of curvature is chosen to 
be 10 inches (k . i/r . q.I) and the torsion t is zero since the beam 
lies in a plane. The beam is uniformly twisted and has a total twist 
of IT radians (in other words, as the beam is traveled from one end to 
the other the cross-section rotates about the elastic axis a total of 

180 degrees). Thus, from equations (5) the curvature components of the 
undeformed beam are 

= (0.1 1n‘^) sin Y 

ky = (0.1 1n‘^) cos Y 

= 0.1 rad/in 

where the angle y defines the orientation of e relative to the 
curve’s normal, n (see figure 2). As shown in figure 5 the beam is 
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assumed to have a rectangular cross-section of dimensions 1 inch by 2 

7 2 

inches. The elastic moduli are taken to be E = 10 Ib/in and 

fi 2 

G = 4 X 10 Ib/in^. The Saint-Venant warping function and torsional 
stiffness constant for a rectangular cross-section as given by Wang in 
reference 27 were used in the computations. 

The solutions were obtained using a computer program based on the 
previously described transfer matrix technique. Twenty-one equally 
spaced stations (including the end-points) along the beam's length were 
utilized. The first four natural vibration frequencies were found to 
be the following: 


Mode 

Frequency 

First 

939 rad/sec 

Second 

3169 rad/sec 

Third 

4490 rad/sec 

Fourth 

7535 rad/sec ^ 


Because of the twist the modal displacements u and v (dis- 
placements in e and e directions) are coupled. The mode shapes 

..y 

are easier visualized by considering the displacements in the directions 
of the normal n and binormal b of the elastic axis (see figure 2). 
Letting n and c be the displacements of the elastic axis in the n 
and b directions, respectively* then n is the inward radial dis- 
placement and c is the out-of -plane displacement. The first mode 
has predominantly out-of-plane motion. The modal displacement c is. 
given in figure 6 as a function of the beam's axial coordinate s. The 
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s, inches 


Figure 6.- Modal displacement c for first mode of twisted curved 
beam example. 
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other displacements are not shown since they are very small compared to 

C- 

Since there is no experimental or calculated data for a beam such 
as the one considered, the results can only be compared to those of an 
untwisted beam. In reference 28 Den Hartog gives expressions for the 
first out-of -plane frequency and for the first in-plane frequency of 
an untwisted cantilevered ring segment. For an untwisted beam having 
the same dimensions as those of the numerical example and having its 
cross-section oriented the same as the one at the ends of the beam of 
figure 5, the first out-of-plane freqiiency is 1049 rad/sec and the 
first in-plane frequency is 5218 rad/sec according to reference 28. 

For the same untwisted beam the analysis and solution method presented 
herein gives the first out-of-plane frequency as 981 rad/sec and the 
first in-plane frequency as 4670 rad/sec. 

Twisted ring .- The natural vibration characteristics were also 
computed for a twisted beam whose elastic axis forms a complete 
circular ring. This beam has the same cross-sectional dimensions and 
the same curvatures as the previous example. The total twist of the 
ring is 2ir radians. Also, forty-one stations were utilized in the 
computations. The first four natural frequencies were found to be 
the following: 
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Mode 

First 

Second 

Third 

Fourth 


Frequency 
1981 rad/ sec 
2307 rad/sec 
5782 rad/sec 
5927 rad/sec 


Figure 7 gives the modal displacements for the first mode. 

As before, the calculated frequencies may be compared with. those 
of an untwisted ring. For an untwisted ring the out-of -plane and in- 
plane motions are not coupled. Consider an untwisted ring having the 
same dimensions as the numerical example and with the longer sides of 
the cross-section parallel to the plane of the ring. According to 
reference 15 the first two out-of-plane and in-plane frequencies for 
the untwisted ring are 


Mode 

First 

Second 


Out-of-plane In-plane 

frequency frequency 


1485 rad/sec 3044 rad/ sec 

4252 rad/sec 8610 rad/sec 


Discussion .- For an untwisted curved beam lying in a plane such 
that ^ ^ twelve governing differential 

equations uncouple into two sets of equations with each set consisting 
of six coupled equations. The variables associated with each set of 
equations are 
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0 47r 8 tt 12tt 16t\ 20tt 


s, inches 

Figure 7.- Modal displacements for first mode of twisted ring example. 
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. (a): u, w, Ojj, M , F^, 

and 

(b): V, tty, Fy 

The variables of set (a) describe the in-plane bending and extensional 
vibration modes. The out-of-plane bending and torsional vibration 
modes are described by the variables of set (b). 

For such a beam the uncoupling of the twelve first-order equations 
can cause trouble in the numerical solution for the modal functions. 

As noted earlier, after determining a frequency the first step for 
obtaining the associated modal displacements is to set one element (or 
variable) of {Y^} equal to unity and then solve equation (150) for the 
remaining elements. For example, if the last element of {Y^}, which 
is F^, is set equal to unity then equation (150) may be written as 



where U-j-j is an 11 x 11 matrix. The first eleven equations contained 
in the above may be written as 

[Uii]{V = -<V 

which can be solved for {Y^}, the remaining elements of {Y^}, For the 
case of the untwisted beam the variable which is set equal to unity 
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must be non-zero. For example, if the frequency has been determined 
for an out-of-plane mode, then the variable which is chosen to be 
equal to unity must be one of the variables of set (b) (setting = 1 

for this case yields modal displacements which do not satisfy the 
boundary conditions). However, the computer program cannot be told 
ahead of time which frequencies will be in-plane mode frequencies and 
which will be out-of-plane mode frequencies. Instead of writing 
separate computer programs for the in-plane and out-of-plane modes, it 
is thought that for cases in which the equations uncouple the problem 
is best remedied by obtaining two modal solutions for each frequency. 
This is done by setting equal to unity, in turn, one variable from 
each of the variable sets (a) and (b). Of course, only one of the 
modal solutions obtained for each frequency is valid. 



IX. SUMMARY 


Developments have been presented for the .extensional equations 
which govern the deformation of curved and twisted space beams. Also, 
discussions of the equations in view of previous work were given. 

First, the exact curvature equations for the deformed beams were 
derived. It is shown that allowing extension of the elastic axis . 
Instead of assuming inextenslon results in an additional term in each 
of the linearized curvature expressions. These equations were used in 
the development of the normal and shearing components of the strain 
tensor. The stress resultants were derived from the strain-displacement 
relations instead of assuming the resultants to be proportional to 
changes in the curvatures. It was shown that the classical approach 
could lead to considerable error in the torsional stress resultant. * 
Also, for the same cross-seetipn the presented stress resultants are 
applicable to beams having larger curvatures than those allowed by the 
classical theory. Moreover, by. developing the stress resultants from 
the strain - displacement relations, additional terms appear due to 
the Initial twist of the beam. A vector derivation of the equilibrium 
equations was given. 

The governing equations for natural vibrations were summarized in 
the form of twelve first - order differential equations. A transfer 

V ■ • 

matrix method was described for obtaining the solutions. Numerical 
examples were presented which illustrate the effect that twist has on 
the natural vibration frequencies of curved beams. 
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